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We calculate the Chern-Simons diffusion rate in a strongly coupled A/" = 4 SUSY Yang-Mills 
plasma in the presence of a constant external U{1)r magnetic flux via the holographic correspon- 
dence. Due to the strong interactions between the charged fields and non-Abelian gauge fields, the 
external Abelian magnetic field affects the thermal Yang-Mills dynamics and increases the diffusion 
rate, regardless of its strength. We obtain the analytic results for the Chern-Simons diffusion rate 
both in the weak and strong magnetic field limits. In the latter limit, we show that the diffusion 
rate scales as B xT^ and this can be understood as a result of a dynamical dimensional reduction. 



I. INTRODUCTION 



The non-Abelian gauge fields possess a rich topological structure. The vacuum contains an infinite number of 
energy-degenerate sectors characterized by an integer: the Chern-Simons number Ncs that is a topological quantity 
determined by the global structure of the gauge fields. At zero temperature, the different Chern-Simons sectors 
are connected by quantum tunneling transitions (the instantons [DH]); this leads to the picture of "9- vacuum" of 
non-Abelian gauge theories. At finite temperature, the gauge configurations that change the Chern-Simons number 
can also be activated thermally. We will, loosely speaking, refer to these configurations as "sphalerons" both at 
weak and at strong coupling, even though this term is usually reserved for describing the classical solutions at weak 
coupling [3113]. As opposed to the tunneling processes (instantons), the sphaleron rate is not necessarily exponentially 
suppressed [311]. 

The change of Chern-Simons number in such a process is given in terms of the topological Pontryagin invariant: 

^Ncs -^Jd'x F^^^Pn^) = £^l d'xtTE ■ B, (1) 

where E and B are non-Abelian electric and magnetic fields, and g is the Yang-Mills coupling. The rate of change of 
Chern-Simons number is called the Chern-Simons diffusion rate Tcs- It is simply the probability of a Chern-Simons 
number changing process to occur per unit volume and per unit time: 

Tcs = = / d^x ( J^F;,Fr(x) J^f^„",Fr^(0)) (2) 

The weak coupling result for the diffusion rate for SU{2) gauge theory is given by [7]: 

Tcs = KV°log(l/.g2)r'* (SU(2), weak coupling) , (3) 

where n' is a (numerically large) constant. This expression can be understood as a result of the Langevin~type 
dynamics of the non-Abelian gauge fields [7] [S] . Since the sphaleron at the peak of the barrier separating Chern- 
Simons sectors is a purely magnetic field configuration, the factor (g^T)^ in this expression can be understood as the 
inverse magnetic screening length that determines the characteristic inverse volume of the sphaleron, and g'^T\og{l/g) 
- as the typical inverse time scale of the process. Numerical calculations show that the constant k' ~ 10 [8l |9] for 
SU{2). 

The holographic AdS/CFT correspondence [T0Hl2] makes it possible to compute the Chern-Simons diffusion rate 
in A/" = 4 SYM plasma at — > oo in the strong coupling regime [13] . The strong coupling result is given by [T3] : 

Tcs = ^LHr^^^ (AA - 4 SYM, strong coupling, large N) (4) 
2567r'^ 

It is remarkable that compared to the weak coupling case, at strong coupling the diffusion rate is substantially 
enhanced. This means that Chern-Simons diffusion is not exclusively the property of semi-classical fields, as the 
(Minkowski boundary) dynamics at strong coupling has to be driven by quantum effects. 
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The Atiyah-Singer index theorem relates the change in Ncs to the chirahty change in the fermionic sector - 
every Chern-Simons number changing transition is accompanied by the flip of chirality. In thermal equilibrium, the 
sphalerons thus lead to the decay of any excess chiral charge = (Jg) present in the medium. A linear response 
relation for a small chiral chemical potential leads to [HI [S] the decay: 

— ^-CN,— (5) 

where the constant C depends on the details of the fermionic sector. 

These TV^s-changing transitions violate locally P and CP symmetries and lead to very interesting consequences in 
a number of various physical settings. In weak interactions, they are related to the violation of baryon plus lepton 
{B + L) number^ and play a crucial role in electroweak baryogenesis scenarios [5| [6t iMlIlQ) . In strong interactions, 
there exists a very strong constraint on the amount of global P and CP violation; this constraint originates mainly from 
the experimental upper bound for the neutron dipole moment |20j . However on theoretical grounds, the observed CP 
invariance is still lacking a conclusive explanation ("the strong CP problem"). On the other hand, the instantons and 
sphalerons provide explicit examples of the local fluctuations of topological charge leading, through the index theorem, 
to the local imbalance of chirality in QCD plasma. In a medium with high enough temperature, the fluctuations might 
lead to observable effects such as the "chiral magnetic effect" [51] which is an induction of an electric current by an 
external magnetic field in the presence of Chern-Simons number changing processes. 

In this paper we address the effect that an external Abelian magnetic field has on the topological fluctuations in 
the plasma. This question has been addressed in the past in several different contexts. The chiral magnetic effect is 
seen in various lattice simulations [22H24] . The effects of magnetic field on an instanton configuration is also studied 
on the lattice [551 HB] • In [13 the Dirac spectrum of an instanton in the presence of magnetic field is analyzed and it 
is argued that in strong magnetic field regime some of the physical quantities such as the magnetic and electric dipole 
moments are dominated by the (near-) zero modes. In [28', it is shown that the diffusion rate due to an electoweak 
sphaleron in the Higgs phase increases with a presence of an external magnetic field. 

Here we calculate the diffusion rate for strongly coupled A/" = 4 plasma in the presence of an external magnetic 
field via holography. The dual gravity is characterized by a full solution of asymptotically AdS^, five-dimensional 
Einstein-Maxwell system with a constant magnetic flux. Since the solution takes into account all of the back-reaction 
of the magnetic held on the metric, this gravity description allows one to work with arbitrarily strong magnetic fields. 
This configuration was previously studied in [29^ and, together with its various extensions, was pursued in a different 
context in 'HDI. 



The rest of the paper is organized as follows. In section 
in [2S1 which constitutes our dual metric. Then, in section 



Hi we summarize the magnetic brane solution presented 



HI we calculate the diffusion rate and present the exact 



numerical result. Section IV is devoted to the high temperature and low temperature limits of this result, for which 
we obtain the analytic expressions. At high temperature we calculate the leading term to the zero magnetic field result 
Q analytically by treating the magnetic field perturbatively. At low temperature, we use the (2-|-l)-dimensional BTZ 
black hole solution to describe the dimensionally reduced metric and express the rate also analytically at this limit. 
We conclude the paper by arguing that the effect is negligible for heavy ion collisions and discuss a possible implication 
for a particular scenario of electroweak baryogenesis. 



II. THE DUAL GEOMETRY 



Let us briefiy review the magnetic brane solution studied in |29j that constitutes the gravity dual of our problem. The 
basic setup is the five-dimensional Einstein-Maxwell theory with a negative cosmological constant. In the boundary 
gauge theory, the Abelian Maxwell field is associated with the U{\)r symmetry which the gauginos, the chiral 
multiplet fermions and scalars are charged under. The action also contains the five-dimensional Chern-Simons term 
that accounts for the U{l)ii anomaly: 

S = -^J d^^V^, [r + ^--Fm. - ^) + ^ JaaFAF + (6) 



We will refer to the violation of {B + L) simply as baryon number violation from now on. 
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The last term in ^ is the boundary term that is determined by imposing a sensible variational principle. The AdS 
radius I will be set to unity for the rest of the paper. It is also worth to mention that the coefficient of the variation 
of the Chern-Simons term is fixed by the U{l)fi anomaly of the gauge theory. This relates the field strength F in ^ 
to the physical field strength in the gauge theory as = ^/SF. We refer the reader to l29] for further details. 

Let us now introduce a constant magnetic flux in the X3 direction: 

F = B dx^ A dx^ (7) 

The constant flux satisfles the Maxwell equations trivially. For the Einstein equations 

Rmn — ^gMN + -^F^^ FABguN — 2FmaF]^ (8) 



we start with the general form of the 5-dimensional metric 

dr_^ 

W) 



ds" = ~U{r)dt' + -^+ e^^'^-\dx\ + dxl) + e^^^^^Uxl (9) 



and seek for asymptotically AdS^ solutions with a horizon to study the problem at nonzero temperature, in accordance 
with the AdS/CFT correspondence. Unfortunately, there is no known analytical solution to the Einstein equations 
satisfying these properties. However it is possible to solve them numerically. Notice that the metric components 
U, V and W are functions of a single variable r. Therefore the Einstein equations ^ are a set of coupled ordinary 
differential equations in r. Following [53], we rescale the coordinates such that the horizon is at r = 1 (i.e C/(l) = 0), 
and [/'(!) = 1,1^(1) = W{1) — 0. The remaining two initial conditions V'{1) and H^'(l) are encoded in the Einstein 
equations. Then, given the initial conditions, we integrate the Einstein equations from the horizon r = 1 to the 
boundary r — > 00 numerically. The result is that the solutions indeed are asymptotically AdS^: 

U{r)^r^, e2^('■)^^;r^ e^^^''^ ^ wr'^ as r ^ 00 (10) 

Here v and w are the functions of the magnetic field strength B. To get the correct AdS 5 limit on the boundary, we 
should rescale the xi,X2 and X3 coordinates: 

^j,2 2V{r) J2W{r) 

ds^ = ~Uir)dt^ + —— + (dxl + dxl) + dxl (11) 

This rescaling also affects the form of the physical magnetic flux: 

F = —dxi A dx2 (12) 

V 

It is also useful to define the physical magnetic field strength in the gauge theory B: 



B = V3- (13) 

V 



III. THE CHERN-SIMONS DIFFUSION RATE 



The diffusion rate ([2| is nothing but the zero frequency and wavelength limit of the symmetrized Wightman 
correlator of the topological charge: 

G^y"^{oj,k) ^ \ J d'xe-^'^^(^l^\F;^^Fr{x) , ^F^.i^rw}) (14) 

The fluctuation-dissipation theorem relates the symmetrized Wightman correlator to the imaginary part of the Green's 
function: 

G'y^'ioj.k) = -coth (-^) Im[G"{iJ,k)] w - — Im[G"{uj,k)] (15) 



where the retarded Green's function is: 

/■ . ., /fl - 1 - 1\ 

(16) 



G^(w,^) = d^xe-''''=e{t) ^ 



\F^Jrix),\F^.Fr{0) 



(17) 
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Hence, to extract the diffusion rate all we need to know is the zero momentum and small frequency limit of the 
retarded Green's function: 



cs 



I lim — Im 

8n^ J w-i-o uj 



G^(a;,fc = 0) (18) 



The imaginary part of the retarded Green's function is associated with the thermal dissipation and therefore is a T- 
odd function. As a result the leading term in the small frequency expansion has to be linear in w, with the coefficient 



that controls the diffusion rate ( 18 ) . 



We compute the retarded Green's function of the topological charge by using the AdS/CFT duality where the 
gravity side is described by the metric ( |ll| ). The duality relates the operator ^F^j^F^^" in A/" = 4 SYM with the 
RR scalar Co which is the axion in the ten-dimensional type IIB supergravity theory. In the absence of the U{l)fi 
Maxwell field, the axion is nothing but a free massless scalar in the consistent truncation to the five-dimensional 
effective theory after Kaluza-Klein reduction of the . To the second order in fluctuations in the fields, it remains 
as a free scalar even with the presence of a background J7(1)_r flux. To see this let us analyze the axion part of the 
ten-dimensional type IIB SUGRA action in the string frame: 

-CoAdSp), (19) 

where Cq is the axion, C2 and B are anti-symmetric RR and NS-NS 2 forms respectively. In the gauge theory, 
the U{1)r symmetry is realized as the diagonal U{l)n of the Cartan subgroup U{l)j^ of SU{4:)]i. In the SUGRA 
picture, this corresponds to wrapping around the 3 angles of the simultaneously by a special choice of truncation 
parameters. The Kaluza-Klein reduction can be written as [511 15^ : 



i=l 



where gMN is the five-dimensional AdS metric, /i^'s are 2 independent the parameters with the constraint X]i=i A^f — 
1 and ^i's are the remaining 3 angles of S^. The Kaluza-Klein reduction of this metric is shown [211 132] to be a 
consistent truncation that generates the action To see that the Maxwell sector does not affect the axion sector 



it is sufficient to observe that the only place where the U{l)ii field might show up in the action (19 1 is the ten- 
dimensional metric determinant \/—G. However this is not possible because the gauge field can be thought of as a 
re-parametrization of the angles (jii of under which G is invariant. From the gauge theory perspective, there cannot 
be any gauge dependent term in the action; therefore G has to be independent of which is in fact true. Also the 
NS-NS 2 form B can be turned off consistently in the compactifiaction. As a result, the axion field Cq is a massless 
scalar in the five-dimensional Einstein-Maxwell theory. 

The associated Klein-Gordon equation for a generic metric of the form ^ is not analytically soluble, since we do 
not even know the metric explicitly. However it can be expressed in terms of the functions C/, V and W that determine 
the metric at the linear order in frequency, which is the limit we are interested in. The dissipative processes that lead 
to an imaginary part in the retarded Green's function are associated with the in-falling boundary conditions for the 
dual scalar field near the black brane horizon. We shall study the Klein-Gordon equation in Eddington-Finkclstein 
coordinates where the in-falling boundary conditions near the horizon appear naturally. Let us define 

r r' 

U=t+ — (21) 
J 00 U{r') 

the metric ^ with this coordinate transformation becomes: 

ds^ = ~U{r)dtl + 2dr dU + e^^^'^^ {dxl + dxj) + e'^^^^'Uxj (22) 

Note that the frequency remains unchanged: lo^ ^ idt, = idt ^ uj- Furthermore, the wavefunctions are of the 
form p-^'^^-'+i-^-^ (f)^ fe('')- The mode function (f)^ ^{r) should have regular near horizon behavior. Near the horizon, 
fa t + ijr^j.^-^ In I?" — '''/il up to constant, where r — is the horizon. The Hawking temperature for a general U{r) is 
given by the relation: U'(rh) — AttT. Therefore the wavefunction has the expected infalling behavior e~''^*|r—r'/i| 
which arises naturally in the Eddington-Finkclstein coordinates. 



^ 47rT 



5 



The Klein-Gordon equation d^{y/—g g^'^d^cj)) = with the metric (22 1 reads: 



C,^+(/W+lnC/(r))'<^:,-^ 



(23) 



Here we defined e^^^'' 



J—g = ^v(r)^'W(r') g^j^^ / jg short for dr- We now take fc = and solve (23) for small 
frequencies. To do this, we expand 4>Lj{r) to linear order in w for u! << 1^. 



^u,k=o = Mr) + iojMr) + 0{^') 
The Klein-Gordon equation is now expanded in the same fashion as: 

0^' + (/ + In - icj -I- ^4,0 - (P'l - (/ + ln;7)V'i 

which should be solved order by order in w. 
The 0{ijJ°) order equation has the solution: 



(24) 



(25) 



4'oi'r) = ci + C2 



dr' 



(26) 



The second term is divergent near the horizon since U{rh) = 0. In order (j)Q{r) to be well behaved near the horizon, 
C2 must vanish. We also set ci=l. This choice translates into normalizing the homogeneous axion source in the 
boundary theory to unity in the Ad S/CFT dictionary. 

After plugging 0o(r) = 1 into (251, the C'(w^) equation becomes: 



(/>'/ + (/ + In [/)Vi 



which can be solved as: 



dr 



J(rh)-f(r')\ _ 



u 



dr' 
U{r') 



1 - 



\/-9{rh) 
v/=5F), 



(27) 



(28) 



Here the constant \/ —g{rh) is chosen to assure the regularity of the solution near the horizon. The other integration 
constant is determined by setting the limit of the integral such that (pi vanishes near the boundary and does not effect 
the normalization of the source. 



As a result the solution of the Klein-Gordon equation (231 in the zero momentum, small frequency limit is: 



\k=0 



= 1 + iuj 



= 1 



dr' ( _ ^/-g{rh) 
U{r') \ y=ffp)^ 

dr 



(29) 



We now expand (29 1 near the boundary. Recall that our metric converges to the AdS^ metric r ^dr^ -I- r^dx^dx^ 
near the boundary r — > oo. Therefore yj—g — > r^, U{r) — > r^ and 



\k=0 



, . iu:^-g{rh) 



(30) 



which reproduces the expected boundary behavior of a massless scalar. The first term is the non-normalizable mode 
that corresponds to the homogeneous axion source 5$ normalized to unity. The coefficient of the term is: 
F^"^), that is the response of the system to the axion source (5$. Putting back all the constants it is found as: 



iuj^J-g{rh) 



(31) 



Here there is a slight abuse of notation. By writing << 1 we actually mean << 1. But it is always possible to work with 

dimensionless frequency by rescaling t» and r which we are doing by setting ?7'(1) = 1 (i.e. 47rT=l). 
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Here A^^, = 4 is the scaling dimension of the topological charge operator. We can easily calculate the retarded 
Green's function in the zero momentum, small frequency limit through the linear response relation (Kubo's formula): 



5l -Ff„F^'''(u 



Thus we obtain in the leading hncar frequency hmit: 



G^(w,fc = 0) = 



(5$ 



Using the properly scaled metric (11) it reduces to: 



(32) 



(33) 



G^(w,fc = 0) = 



lUJK 



Vy/W 

Recall that the diffusion rate is related to the retarded Green's function as: 

.2x2 



(34) 



res = 



9_ 
87r2 



g 

87r2 



lim — Xm 
2T 



G^{u,k = 0) 



37r^ vww 



(35) 



In the last step we have used the holographic relation between the gravitational and gauge couplings: k 

— i^ where Gr is the five-dimensional Einstein constant. Recall that the functions v and w are obtained in the 
units C/(l) = 0, U'{1) = 1. In these units T = l/iir and due to conformal symmetry we can choose our physical 
parameter to be the dimensionless magnetic field strength B/T^. We then express the diffusion rate in terms of this 
dimensionless ratio as follows: 



r(S/T2 



TiB,T) 



To 



To = 



2567r3 



(36) 



This quantity measures the ratio of the diffusion rate in a medium with a magnetic flux to the one in a medium 
without a magnetic flux, Tq. From the plot (Fig ij, we see that T{B/T'^) is a monotonously increasing function of 
the magnetic field strength. We now investigate the high and low temperature limits of the diffusion rate. 



IV. HIGH AND LOW TEMPERATURE LIMITS 



A. Weak magnetic field (high T) limit 

In the limit where B << T^, the magnetic field can be treated as a perturbation to the T"* scaling of the diffusion 
rate. Due to the R-charge symmetry of the gauge theory, we expect the perturbation series to start at the order B^. 
Also we can argue that we do not expect the result to depend on the direction of the magnetic field at the leading 
order and the analyticity of the perturbation series dictates the B^ dependence. It is also clear from the gravity 
side that the perturbative corrections to the Einstein equations (|8| will be at the order B^ since the magnetic field 
enters the equations in a gauge invariant way. Let us perturb the Einstein equations ([s]) around the AdS^ black brane 
solution: 

U{r) = r^-'^+B^uiir) 

V{r) ^\og{r) + B^ vi{r) 
Wir) = log(r) + B^ Wi (r) (37) 

where th = ttT. To extract the diffusion rate, wc only need to know y/^^g at the perturbed horizon. Let us call the 
perturbed horizon rp, such that U{rp) = 0. By expanding to first order in B^ we get: 

= e2yM+w(r) =4(1+ B\2v,{rp) + wi{rp)) (38) 
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f(S/T'^2) 




FIG. 1: The Chern-Simons diffusion rate T = r(Z3, r)/ro(T), normalized by the zero magnetic flux value To = 2^^T^ as 
a function of the dimensionless magnetic field B/T'^ . The diffusion rate is monotonously increasing with the magnetic field 
strength and has the asymptotic behavior F ~ BT^ for B >> . 



From the linearized (in B ) Einstein equations we obtain the following relations: 



2 
972 

2v[{r) + w[{r) ^ (39) 



These corrections should vanish as r — > 00 since we should have asymptotically AdS^ behavior. This condition sets 
the constant 2vi{r) + Wi(r) — leading to: 

+ + ,40) 



As a result, to leading order in we get 

2567r3 " V" ' Gn* ' "^T^ 



Comparing with the exact numerical result ( 36 1 we have found that this approximation differs from the exact result 
only by 2.5% for B ^ 10 . This suggests that a proper expansion parameter is probably B'^/{ttT)*. 



B. Strong magnetic field (low T) limit and BTZ black hole 

For strong magnetic field B >> , v and w have the limits 

lim = — ,^ = ^ (42) 

B^oo 6w(47rr)2 96\/37r2T2 



fn the last step we plugged in the physical magnetic field strength B — ^/iB/v. Therefore in the presence of a strong 
magnetic field the diffusion rate is: 

T{B,T)^^^^^BT^ ' (43) 

It is possible to obtain the same result in a way where the physics is more transparent. In |29) . it was shown that 
there is another solution to Einstein-Maxwell configuration which is a product of a (2-|-l)-dimensional BTZ black hole 
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in the (r, t, X3) directions and a fiat surface in the transverse plane: 



ds^ = -3(r2 - rl)dt^ + f' + ir^dxl + ^{dxl + dxl) (44) 

The Hawking temperature is T = 3r/i/27r. This is an appropriate description of the dual gauge theory in the strong 
magnetic field limit. When ^/B » T, the fields occupy the lowest Landau levels and transitions to higher Landau 
orbits are suppressed by the factor e^^^^ . As a result the motion in the transverse plane is frozen and the only 
remaining degree of freedom is along the longitudinal direction X3 where the motion is not affected by the magnetic 
field. The overall factor of B in front of the transverse plane corresponds to the density of the lowest Landau levels in 



the gauge theory. The geometry ( 44 ) is therefore merely the holographic manifestation of the dimensional reduction 



of the system in the presence of a strong magnetic field. Using -\/— 5 = r^B and the expression (33) for the retarded 
Green's function we obtain 

G^(a;, fc = 0) = iujKBrh = iujK — + 0{lj'^). (45) 

3v 3 



Note that the expression (33) is valid only for metrics that are asymptotically AdS^, since the retarded Green's 
function is read off from the coefficient of r^^ in the expansion of 4> near the boundary. Here in the dimensionally 
reduced calculation, we assume that the magnetic field strength is very close to UV scale of the gauge theory so that 
the transition from BTZ to AdS^ occurs close to the UV boundary and far away from the horizon in the dual picture. 
Therefore, the BTZ approximation is good near the horizon and we can use the BT Z m etric for \/—g{rh)- 



Plugging this result for into the expression for the diffusion rate (first line of (36)) leads to the correct strong 



magnetic field result (43). This result can be interpreted as follows: Tcs is a transport coefficient and its dynamics 
are governed by the IR sector of the theory. In the R -charged sector, the IR dynamics is dominated by the lowest 
Landau level fermions. The contribution of scalars to IR is suppressed exponentially since they have no zero modes 
as opposed to fermions. The holographic dimensional reduction mechanism explained above suggests that the strong 
interactions between the dimensionally reduced R-charged sector and R-neutral non-Abelian gauge sector effectively 
induces a dimensional reduction for the gauge field in the IR, even though it does not see the magnetic field at the 
level of Lagrangian. Therefore the Ghern-Simons number changing transitions occur in one spatial dimension which 
brings a factor of to the rate on dimensional grounds. This should be multiplied by the density of lowest Landau 
levels in transverse plane which scales as B/tt. 



V. CONCLUSION AND DISCUSSIONS 



In the strongly coupled Af — 4: theory, the existence of an external U{1) u magnetic field increases the Ghern-Simons 
diffusion rate in comparison to the zero magnetic field case. This modification is due to the strong interactions between 
the R-charged sector of the theory and the R-neutral Yang-Mills fields that is realized through the alteration of the dual 
gravity by the magnetic field. Furthermore in the strong magnetic field limit the Ghern-Simons number transitions 
occur in the dimensionally reduced system and the rate thus scales as B xT"^. The absence of the interaction between 
the C/(l)/j and Yang-Mills fields makes the ?7(1)_r magnetic field a reasonable proxy for a "real" U(l) magnetic field. 
Wc thus expect a similar behavior in a strongly coupled, non- SUSY gauge theory with a U{1) magnetic field. 

We would like to point out an interesting analogy between the electroweak sphaleron in a magnetic field and our 
result. The electroweak sphaleron in the Higgs phase (Klinkhamer-Manton sphaleron) [3] is spherically symmetric 
for zero Weinberg angle. However a non-zero Weinberg angle reduces the symmetry into an axial one [331 I34j . 
Furthermore, the deformed sphaleron develops a magnetic moment. Due to this dipole moment, an external magnetic 
field decreases the energy of the sphaleron which lowers the barrier between the topologically inequivalent vacua 
PSI [55] . This puts some constraints on the electroweak baryogenesis scenario with a primordial magnetic field since a 
lower barrier means topological transitions are more frequent and they wash out any initial baryon number asymmetry. 

Our result suggests that an external magnetic field raises the rate of topological transitions also in strongly coupled 
regime. In this case the energy of the sphaleron is altered by the magnetic field due to strong interactions. Furthermore, 
the spherical symmetry is reduced to axial symmetry through dimensional reduction since there is a residual 0(2) 
symmetry in the transverse plane. These effects of the magnetic field are, of course, due to a different mechanism 
than the Kinkhamer- Manton sphaleron case above. However, it is noteworthy that even though the mechanisms are 
quite different, it seems a coupling to an extra C/(l) field affects the non-Abelian sector in a similar fashion for both 
cases. 
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There are electroweak baryogenesis scenarios where there exist CP violating processes in the scalar (Higgs) sector 
[I4l - [T7] . These processes can be described by an effective inhomogeneous axion field where chiral quarks are created 
locally on these inhomogeneities. These inhomogeneities correspond to the walls of the bubbles of broken phase 
nucleated inside the symmetric phase. As the created net chirality propagates into the symmetric phase it triggers the 
electroweak sphaleron transitions in the symmetric phase which ultimately creates the baryon number. However, the 
strong sphaleron transitions damp the net chirality and therefore reduce the baryon number creation [9l [18] . Taking 
into account the existence of a primordial magnetic field, basing on our result we can conclude that it would further 
decrease the rate of baryon number generation if the relevant strong interactions are in the strongly coupled regime. 
Therefore, it seems that an existence of a primordial magnetic field works against the baryon number generation. 

Finally, let us pass from the big bang to the little bang: the heavy ion collisions. The background magnetic field 
created in these collisions is typically of the order T ~ ^/B ~ [HJ [35] . Our result suggests that a magnetic field 
of this strength has a negligible effect on the topological charge diffusion. To be more precise 

T{B,0) 67r4 ^ ' 

for B — T^. Therefore the magnetic field changes the rate only by 0.17 %. This means that the quantitative estimates 
of the chiral magnetic effect in heavy ion collisions can be safely done by using the sphaleron rate in the absence of 
magnetic flux. 
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